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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 


U.S.  customary  units 
metric  (SI)  units  as 

of  measurement 
follows: 

used  in 

this  report  can  be  converted  to 

Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  Inches 

•  6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x 

10'3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

lTo  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  formula:  C  *  (5/9)  (F  -32). 

To  obtain  Kelvin  (K)  readings,  use  formula:  K  «  (5/9)  (F  -32)  +  273.15. 
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SYMBOLS  AND  DEFINITIONS 


ajm/k  roughness  term 

Cg  group  velocity  (meters  per  second) 

D  damage  to  the  profile  (cubic  centimeters) 

D'  dimensionless  damage 

d  water  depth  in  flat  part  of  tank  (centimeters) 

d 1 5  equivalent  diameter  of  the  riprap  stone;  15  percent  of  the  total  weight 
of  the  armor  gradation  is  contributed  by  stones  of  lesser  weight 
(millimeters) 

d05  equivalent  diameter  of  the  filter  stone;  85  percent  of  the  total  weight 
of  the  filter  gradation  is  contributed  by  stones  of  lesser  weight 
(millimeters) 

g  acceleration  due  to  gravity  (9.81  meters  per  second) 

H  incident  wave  height  (centimeters) 

wave  height  at  the  zero-damage  level  (centimeters) 

L  wavelength  (centimeters) 

Lq  deepwater  wavelength  (feet) 

Ng  stability  number 

Nz  stability  number  at  the  zero-damage  level 

NZp  stability  number  at  the  zero-damage  level  divided  by  the  average  of  the 

prototype  stability  numbers  at  the  zero-damage  level 

Rg  Reynolds  number  with  roughness  term 

Rjj  Reynolds  number  using  wave  height 

Rz  runup  at  the  zero-damage  level 

S  distance  from  wave  blade  to  toe  of  embankment  (meters) 

T  wave  period  (seconds)  j 

Wjj  wave  burst  duration  (seconds) 

W50  median  armor  stone  weight  (kilograms);  weight  of  stone  where  50  percent 
of  the  total  weight  of  the  armor  gradation  is  contributed  by  stones  of 
lesser  weight 
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SYMBOLS  AND  DEFINITIONS — Continued 

unit  weight  of  riprap;  2707.1  kilograms  per  cubic  meter 
in  this  study 

unit  weight  of  water;  1000  kilograms  per  cubic  meter 
surf  parameter  (centimeters  squared  per  second) 
angle  formed  between  embankment  slope  and  horizontal 
kinematic  viscosity 
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RIPRAP  STABILITY  SCALL  EFFLCTS 


bli 

Laurie  L.  brodevick  and  John  l.  Ahrens 

I.  INTRODUCTION 

Small-scale  wave  tank  tests  of  riprap  stability  were  conducted  at  the  U.S. 
Army  Coastal  Engineering  Research  Center  (CERC),  and  the  results  were  compared 
with  previously  conducted  large-scale  tests  of  riprap  stability  (Ahrens,  1975) 
to  determine  the  nature  and  magnitude  of  scale  effects.  The  large-scale  tests, 
conducted  in  CERC's  large  wave  tank  (LWT),  used  wave  heights  which  exceeded 
1.5  meters  in  some  instances  and  can  be  regarded  as  prototype  scale.  The  small 
scale  tests  replicated  the  LWT  tests  at  a  1:10  (model:prototype)  Froude  scale. 
Both  small-scale  and  LWT  tests  were  conducted  using  monochromatic  waves  (waves 
of  constant  height  and  period).  The  results  of  the  experiments  were  used  to 
evaluate  scale  effects  correction  factors.  These  results  will  also  be  used  to 
evaluate  future  small-scale  tests  of  riprap  stability  using  irregular  waves, 
the  next  phase  in  this  study. 

It  has  only  been  in  the  last  few  years  that  irregular  wave  conditions  could 
be  satisfactorily  generated  in  the  laboratory  but  only  at  small  scales,  because 
of  the  unavailability  of  a  prototype-scale  irregular  wave  research  facility. 
This  investigation  of  scale  effects  will  allow,  with  some  confidence,  the 
extrapolation  of  the  small-scale  test  results  to  prototype  scale  when  the 
model/prototype-scale  ratio  is  1:10.  It  will  also  give  general  insight  into 
the  nature  of  scale  effects  to  be  expected  when  conducting  model  experiments  at 
other  scale  ratios. 


II.  TEST  SETUP  AND  PROCEDURE 
1.  Large  Wave  Tank  (LWT)  Tests. 

Ahrens'  (1975)  large-scale  tests  were  conducted  in  the  LWT  which  is  193.6 
meters  long,  4.6  meters  wide,  and  6.10  meters  deep.  A  Stillwater  depth  of 
4.6  meters  was  used  for  all  tests.  The  distance  between  the  toe  of  the  embank¬ 
ment  and  the  mean  position  of  the  wave  generator  blade  varied  from  about  119  to 
137  meters,  depending  on  the  slope  of  the  embankment  being  tested.  Details  on 
the  wave  tank  and  generator  are  given  in  Coastal  Engineering  Research  Center 
(1980). 

The  embankment  was  made  up  of  core  material,  a  filter  layer,  and  an  armor 
layer.  The  core  material  was  compacted  bank-run  gravel,  graded  to  the  desired 
slope,  and  was  essentially  impermeable  to  wave  penetration.  The  filter  layer, 
15  to  21  centimeters  thick,  was  placed  between  the  core  material  and  the  armor 
layer.  The  filter  stone  was  sized  such  that  the  ratio  of  the  15  percent  finer 
diameter  of  the  riprap  stone,  djs,  to  the  85  percent  finer  diameter  of  the 
filter  stone,  dgs,  was  usually  less  than  4  and  always  less  than  5.  The  armor 
stone,  which  was  a  diorite  with  a  specific  gravity  of  2,71,  was  divided  into 
three  stockpiles  according  to  the  median  weight,  Wsq*  One  stockpile  ranged 
from  12.2  to  16.3  kilograms,  another  from  33.1  to  35.4  kilograms,  and  the  third 
was  constant  at  54.4  kilograms.  The  relative  size  gradations  of  the  three 
stockpiles  were  the  same;  the  specified  maximum  stone  weight  was  four  times  the 
median  weight  and  the  minimum  weight  was  one-eighth  of  the  median  weight. 
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In  the  LWT  tests  the  following  parameters  were  varied  systematically:  wave 
height,  embankment  slope,  riprap  weight,  and  wave  period.  Wave  heights  varied 
from  0.43  to  1.83  meters;  wave  periods  ranged  from  2.8  to  11.3  seconds;  the 
embankment  slopes  tested  were  1  on  2.5,  1  on  3.5,  and  1  on  5;  and  the  median 
riprap  weight  varied  from  12.2  to  54.4  kilograms. 

2.  Small-Scale  Tests. 

The  small-scale  tests  were  run  in  one  of  CERC's  small  wave  tanks,  0.46 
meter  wide  by  0.91  meter  deep  by  45.7  meters  long,  which  was  used  to  repli¬ 
cate  the  LWT  at  a  1:10  Froude  scale.  The  width  of  the  small  tank  is  one-tenth 
of  the  LWT,  and  the  distance  from  the  toe  of  the  embankment  to  the  mean  position 
of  the  blade  was  made  one-tenth  of  the  distance  in  the  LWT.  Figure  1  shows  a 
profile  view  of  the  tank  used  in  the  small-scale  tests. 

In  the  small-scale  tests  the  wave  height  and  wave  period  were  varied,  but 
the  embankment  slope  and  median  riprap  weight  were  fixed.  The  embankment 
slope  for  the  small-scale  tests  was  1  on  3.5  which  was  one  of  the  slopes 
tested  in  the  LWT.  The  median  riprap  weight  was  fixed  at  0.034  kilogram  which 
replicates  the  stockpile  of  riprap  with  a  median  riprap  weight  of  34  kilograms 
used  in  the  LWT  when  scaled  down  using  the  Froude  scale. 

The  riprap  armor  layer  in  the  small-scale  tests  was  the  same  material  and 
from  the  same  quarry  as  that  used  in  the  prototype  tests,  diorite  with  a 
specific  gravity  of  2.71.  The  gradation  of  the  model  armor  ranged  in  weight 
from  four  times  to  one-eighth  the  median  weight  of  0.034  kilogram,  the  same 
gradation  as  that  used  in  the  prototype  (Fig.  2).  A  gradation  analysis  was  run 
both  before  and  after  testing  to  determine  if  the  gradation  of  the  armor  unit 
changed  over  time.  As  shown  in  Figure  2,  the  two  gradations  appear  to  be  about 
the  same. 

The  filter  layer  consisted  of  small  gravel  with  a  3-  to  8-millimeter  diam¬ 
eter  and  was  constructed  to  approximately  one-tenth  prototype  scale;  model 
and  prototype  filter  layer  gradations  are  shown  in  Figure  3. 

The  core  material  was  compacted  sand  with  a  median  diameter  of  0.2  milli¬ 
meter.  No  attempt  was  made  to  replicate  the  core  material  at  a  1:10  Froude 
scale  which  is  effectively  impermeable  in  both  model  and  prototype. 

Waves  were  run  in  bursts  of  short  duration  with  an  interval  of  about  1 
minute  between  bursts  to  allow  the  wave  energy  in  the  tank  to  dampen  out.  The 
duration  of  the  wave  burst  was  set  equal  to 

Wd  =  ^  (1) 

where  W,  is  the  wave  burst  duration,  S  the  distance  from  wave  blade  to  toe 
of  embankment,  and  Cg  the  group  velocity  of  the  waves.  The  number  of  wave 
bursts  run  at  a  particular  wave  height  in  the  model,  which  was  equal  to  the 
number  of  wave  bursts  in  the  prototype,  was  normally  enough  to  ensure  that  the 
riprap  profile  was  at  equilibrium  for  the  given  wave  height  and  period.  The 
minimum  number  of  waves  run  at  a  particular  wave  height  ranged  from  340  for 
the  longest  period  waves  to  1,050  for  the  shortest  period  waves.  After  the 
required  number  of  waves  had  been  generated,  the  condition  of  riprap  surface 
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was  documented;  the  wave  height  was  then  increased  approximately  10  percent, 
and  wave  bursts  were  generated  again.  This  procedure  was  continued  until  a 
wave  height  was  reached  which  caused  failure  of  the  riprap.  Failure  was 
defined  as  the  riprap  being  shifted  enough  to  expose  part  of  the  filter  layer 
which  was  removed  by  the  wave  action. 

The  incident  wave  height  was  calculated  using  strip-chart  recordings  from 
two  resistance-type  wave  gages  spaced  one-quarter  of  a  wavelength  apart  and 
placed  as  close  to  the  wave  generator  as  feasible  such  that  the  waveform  sta¬ 
bilizes.  The  gages  were  placed  as  close  to  the  blade  as  convenient  to  increase 
the  recording  time  of  the  gages  before  the  waves  were  reflected  from  the 
structure  and  returned  to  the  gages.  The  average  wave  height,  which  was  the 
average  of  the  wave  heights  for  each  gage,  was  measured  by  visual  inspection  of 
the  strip-chart  recordings. 

In  the  LWT  tests,  a  correction  was  applied  to  the  wave  height  because  of 
the  last  wave  effect  (Madsen,  1970).  The  term  "last  wave  effect"  refers  to  the 
occurrence  of  one  to  three  waves  noticeably  higher  than  the  modal  wave  height; 
the  number  of  higher  waves  is  related  to  the  water  depth-to-wavelength  ratio, 
d/L.  This  term  was  used  because  the  highest  wave  usually  occurred  near  the  end 
of  the  burst.  A  high  wave  also  occurred  near  the  beginning  of  the  burst,  causing 
the  highest  waves  to  bracket  the  smaller  waves  of  almost  uniform  height.  These 
smaller  waves  were  considered  the  modal  waves. 

The  highest  waves  in  the  burst  caused  more  stone  movement  than  the  modal 
waves;  thus,  a  correction  was  made  to  the  modal  wave  height.  The  use  of  the 
modal  wave  height  to  characterize  the  height  of  a  wave  burst  would  result  in 
an  invalid  comparison  of  riprap  stability  for  tests  with  different  wave  periods. 
The  correction  to  the  modal  wave  height  was  determined  by  the  depth-to-wavelength 
ratio,  d/L:  a  decrease  from  1.11  for  a  wave  period  of  2.8  seconds  to  1.04  for  a 
wave  period  of  11.3  seconds. 

The  last  wave  effect  was  not  apparent  in  the  small-scale  tests  because  of 
the  initial  and  final  position  of  the  wave  generator  blade.  The  initial  and 
final  position  of  the  blade  in  the  LWT  tests  was  in  the  center  of  its  total 
stroke  where  the  water  particle  velocities  were  at  a  maximum,  creating  some 
irregularities  in  the  initial  and  final  waves.  In  the  small-scale  tests  the 
blade  started  and  stopped  in  a  rear  position  where  the  water  particle  veloci¬ 
ties  were  zero,  causing  no  apparent  irregularities  in  the  wave  burst. 

The  apparatus  used  to  survey  the  filter  layer  and  riprap  armor  layer  con¬ 
sisted  of  six  vertical  sounding  rods  mounted  on  a  rack  that  moved  along  rails 
mounted  on  the  tank  walls.  Attached  to  the  end  of  each  survey  rod  by  a  ball 
and  socket  joint  was  a  foot  measuring  1.8  centimeters  in  diameter  which  was 
approximately  one-tenth  of  that  used  in  the  LWT  tests.  The  model  surface  was 
surveyed  in  the  same  manner  as  that  of  the  prototype.  Surface  elevations  were 
measured  over  square  grid  points  61  by  61  centimeters  (prototype,  one-tenth  of 
that  for  the  model)  apart  on  a  horizontal  plane. 

The  following  procedure  for  the  small-scale  tests  was  the  same  as  for  the 
prototype  tests  (Ahrens,  1975): 

(a)  Place  and  compact  core  material. 


(b)  Place  and  smooth  filter  layer  material. 


(c)  Survey  filter  layer  surface. 

(d)  Place  riprap  armor  stone  by  dumping  from  a  hand-held  can  to 
simulate  the  prototype  procedure  of  dumping  from  a  skip. 

(e)  Survey  riprap  armor  layer  surface  (reference  survey) . 

(f)  During  the  generation  of  the  predetermined  number  of  wave 
bursts,  collect  wave  data  and  visually  observe  the  behavior  of  the 
riprap  and  wave  runup  on  the  riprap  surface. 

(g)  Survey  riprap  armor  layer. 

(h)  Increase  wave  height  approximately  10  percent. 

(i)  Repeat  steps  f,  g,  and  h  until  failure. 

(j)  Conduct  final  riprap  survey. 

III.  METHOD  OF  DATA  ANALYSIS 

Damage  to  the  riprap  armor  layer  was  quantified  by  comparing  the  profile  of 
the  riprap  armor  layer  taken  at  some  wave  height  (damage  profile)  with  the 
profile  taken  before  any  waves  had  attacked  the  riprap  armor  (reference  profile). 
The  comparison  is  shown  schematically  in  Figure  4.  The  change  in  the  reference 
profile  typically  consisted  of  an  erosion  zone  and  an  accretion  zone,  as  shown 
in  Figure  4.  The  volume  per  unit  length  of  the  erosion  zone  was  used  to  quan¬ 
tify  the  extent  of  damage  to  the  riprap,  D.  Using  the  median  stone  weight, 

W50,  to  characterize  the  size  of  the  riprap,  the  dimensionless  damage,  D', 
is  given  by 

D'  =^W50j2/3  (2) 

where  wr  is  the  unit  weight  of  the  riprap  stone;  i.e.,  D'  is  the  equivalent 
number  of  median  size  stones  removed  by  wave  attack  per  median  stone  length. 

The  word  equivalent  is  used  because  D*  includes  about  40  percent  void  spaces. 

The  incident  wave  height  was  made  dimensionless  through  the  use  of  the 
stability  number,  Ns,  which  was  developed  in  Hudson's  (1958)  study  of  the 
stability  of  rubble-mound  breakwaters.  The  stability  number  is  given  by 


where  H  is  the  incident  wave  height,  and  w„  the  weight  of  water.  Since 
freshwater  (ww  -  1000  kilograms  per  cubic  meter)  and  the  density  of  the  stone 
(wr  ■  2707.1  kilograms  per  cubic  meter)  were  the  same  in  both  prototype  and 
small-scale  tests,  (wr/ww  -  1)  ■  1.71  for  all  tests. 

Data  from  one  small-scale  test  (SET-1)  and  one  prototype  test  (SPL-19)  are 
used  in  Figure  5  to  illustrate  typical  damage  trends  observed  in  this  study. 
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Reference  Profile 


Figure  4.  Riprap  damage  profile. 


Figure  5.  Typical  small-scale  and  prototype  damage 
trends  for  d/gT2  ■  0.0144. 
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The  SET-1  test  replicated  the  LWT  SPL-19  test  at  a  1:10  Froude  scale.  Damage 
trend  refers  to  the  Increasing  cumulative  damage  with  Increasing  wave  height. 
Fitted  to  the  data  for  both  tests  are  curves  of  the  form 

D'  -  aNg  (4) 

where  D*  is  the  dimensionless  damage,  Ns  the  dimensionless  wave  height  stabil¬ 
ity  number,  and  a,b  the  dimensionless  regression  coeff icients.  Figure  5  shows 
that  the  regression  curves  fit  the  data  well,  particularly  at  low  levels  of 
damage,  and  provide  a  convenient  method  of  defining  the  damage  trend.  Also 
shown  in  Figure  5  is  the  zero-damage  level  used  in  this  study  (i.e.,  D*  =  2.0). 
D'  =  2.0  is  about  the  lowest  level  of  damage  that  can  be  consistently  detected 
in  the  inherent  scatter  in  the  survey  data. 

The  two  tests  in  Figure  5  had  a  relative  water  depth  of  d/gT2  *  0.0144, 
where  d  is  the  water  depth  in  the  tank,  T  the  wave  period,  and  g  the 
acceleration  of  gravity.  In  comparing  damage  trends,  both  the  prototype  and 
small-scale  tests  were  grouped  by  relative  depth  (see  Table  1)  to  eliminate 
the  possible  influence  of  wave  period  effects  (Ahrens  and  McCartney,  1975). 

Curves  of  the  form  of  equation  (4)  were  fitted  to  the  model  and  prototype 
data,  and  the  following  two  parameters  were  chosen  to  characterize  the  damage 
trend  (Table  1):  the  stability  number,  Nz,  for  D'  =  2.0  which  characterizes 
the  zero-damage  level,  and  the  regression  coefficient  b  (eq.  4)  which  charac¬ 
terizes  the  rate  of  increase  in  damage  with  increasing  wave  height.  The 
parameters  Nz  and  b  are  tabulated  and  grouped  by  relative  depth  in  Table  1 
to  facilitate  comparison  of  small-scale  and  prototype  values  for  similar  wave 
conditions. 


IV.  COMPARISON  OF  MODEL  AND  PROTOTYPE  DATA 

1.  Damage . 

The  values  of  Nz  and  b  from  Table  1  are  plotted  versus  d/gT2  in  Figure 
6  to  demonstrate  the  influence  of  both  the  scale  effects  and  the  wave  period 
effects  on  Nz  and  b.  The  figure  shows  that  the  small-scale  tests  had  lower 
values  of  Nz  and  generally  lower  values  of  b  than  the  prototype  tests  with 
similar  wave  conditions.  This  finding  indicates  that  damage  is  initiated 
earlier  in  the  small-scale  tests  than  in  the  prototype  tests  but  proceeds  at  a 
slower  rate,  with  respect  to  increased  wave  height.  The  convergence  in  the 
damage  trends  typical  of  small-scale  and  prototype  tests  can  be  seen  in  Figure 
5.  The  regression  curves  cross;  however,  the  actual  data  indicate  that  while 
the  damage  levels  in  the  small-scale  tests  may  approach  those  of  the  prototype, 
they  do  not  surpass  them  for  similar  values  of  the  dimensionless  wave  height. 
Figure  7  is  similar  to  Figure  5  except  it  shows  all  the  data  for  tests  where 
d/gT2  =  0.0144,  which  includes  the  data  in  Figure  5.  The  small-scale  and 
prototype  data  fields  overlap  somewhat,  but  the  crossover  suggested  by  the 
regression  curves  in  Figure  5  does  not  occur.  For  the  tests  where  d/gT2  * 
0.0264  and  0.0065,  there  is  more  overlap  or  convergence  of  small-scale  and 
prototype  data  fields  than  shown  in  Figure  7;  for  tests  where  d/gT2  =  0.0037 
there  is  no  overlap  and  little  convergence  in  the  damage  trends.  The  reason  for 
the  convergence  of  damage  trends  is  unclear,  but  it  may  reflect  the  influence 
of  breaker  characteristics  or  may  be  caused  by  the  size  of  the  data  set  and  the 
Inherent  scatter  in  the  data.  Convergence  in  the  damage  trend  indicates  a 
reduction  in  scale  effects  from  the  zero-damage  level. 
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Table  1.  Basic  data 


Test 

designation1 

d/gT2 

— 

T 

U) 

Armor 

layer 

thickness 

(cm) 

Wso 

(ke) 

Nz 

N*P 

b 

«z 

(cm) 

SPL-17 

0.0595 

2.80 

41.76 

34.02 

2.834 

2.648 

6.333 

112.56 

SPL-27 

0.0595 

2.80 

29.56 

12.25 

2.462 

2.648 

5.151 

69.56 

SET-7 

0.0588 

0.89 

5.12 

0.034 

1.870 

— 

3.410 

7.44 

SET-8 

0.0588 

0.89 

4.97 

0.034 

1.644 

— 

2.210 

6.52 

SET-22 

0.0588 

0.89 

4.97 

0.034 

2.267 

6.313 

9.02 

SET-24 

0.0588 

0.89 

4.48 

0.034 

2.604 

5.694 

10.33 

SPL-18 

0.0264 

4.2 

39.32 

34.02 

2.455 

2.136 

4.355 

97.51 

SPL-25 

0.0264 

4.2 

52.43 

34.02 

2.185 

2.136 

5.969 

86.78 

SPL-28 

0.0264 

4.2 

30.48 

12.25 

1.862 

2.136 

4.686 

52.61 

SPL-35 

0.0264 

4.2 

55.47 

54.43 

2.042 

2.136 

3.736 

94.85 

SET-2 

0.0264 

1.33 

5.40 

0.034 

1.816 

3.591 

7.22 

SET-4 

0.0264 

1.33 

4.12 

0.034 

1.772 

3.902 

7.04 

SET-9 

0.0264 

1.33 

5.03 

0.034 

1.610 

— 

3.310 

6.40 

SET-19 

0.0264 

1.33 

5.55 

0.034 

1.431 

— 

2.344 

5.67 

SPL-19 

0.0144 

5.7 

46.33 

34.02 

1.979 

1.970 

7.149 

78.61 

SPL-23 

0.0144 

5.7 

44.81 

34.02 

1.803 

1.970 

5.179 

71.60 

SPL-29 

0.0144 

5.7 

29.87 

12,25 

1.780 

1.970 

6.841 

50.29 

SPL-32 

0.0144 

5.7 

42.37 

54.43 

2.359 

1.970 

8.206 

109.58 

SPL-36 

0.0144 

5.7 

42.06 

54.43 

1.928 

1.970 

6.120 

89.55 

SET-1 

0.0144 

1.80 

4.78 

0.034 

1.490 

— 

3.512 

5.88 

SET-3 

0.0144 

1.80 

4.51 

0.034 

1.412 

4.709 

5.61 

SET-10 

0.0144 

1.80 

5.06 

0.034 

1.605 

5.682 

6.37 

SET-20 

0.0144 

1.80 

4.82 

0.034 

1.624 

— — 

4.083 

6.46 

SPL-20 

0.00646 

8.5 

43.89 

34.02 

2.051 

2.059 

7.614 

81.47 

SPL-24 

0.00646 

8.5 

46.63 

34.02 

2.020 

2.059 

6.611 

80.22 

SPL-30 

0.00646 

8.5 

25.91 

12.25 

2.165 

2.059 

12.195 

61.17 

SPL-33 

0.00646 

8.5 

45.72 

54.43 

2.102 

2.059 

7.573 

97.66 

SPL-37 

0.00646 

8.5 

51.51 

54.43 

1.959 

2.059 

5.396 

91.01 

SET-5 

0.00646 

2.69 

4.45 

0.034 

1.406 

2.422 

5.58 

SET-6 

0.00646 

2.69 

4.36 

0.034 

1.829 

5.040 

7.13 

SET-11 

0.00646 

2.69 

4.48 

0.034 

1.805 

4.856 

7.16 

SET-21 

0.00646 

2.69 

4.54 

0.034 

1.802 

— 

4.705 

7.16 

SPL-22 

0.00365 

11.3 

47.24 

34.02 

2.387 

2.370 

8.145 

94.79 

SPL-31 

0.00365 

11.3 

29.26 

12.25 

2.235 

2.370 

6.687 

63.15 

SPL-34 

0.00365 

11.3 

48.46 

54.43 

2.487 

2.370 

8.480 

115.52 

SET-12 

0.00366 

3.57 

3.84 

0.034 

1.694 

- - 

7.000 

6.64 

SET-13 

0.00366 

3.57 

4.18 

0.034 

1.649 

- - - 

5.525 

6.55 

SET-23 

0.00366 

3.57 

4.82 

0.034 

1.996 

“““ 

8.479 

7.92 

!SPL:  LWT  test;  SET:  small-scale  test. 
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Mean  R/H 


Figure  6.  Zero-damage  stability  numbers,  damage  rate 
coefficients,  and  relative  runup  versus 
relative  depth. 
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Dimensionless  Wove  Height  (Ns) 


Figure  7.  Small-scale  and  prototype  damage  trends 
for  d/gT2  -  0.0144. 


In  addition  to  the  information  on  damage  trends.  Figure  6  shows  that  the 
small-scale  tests  exhibit  less  influence  of  wave  period  at  the  zero-damage  level 
than  the  prototype  tests.  The  prototype  Nz's  have  the  characteristic  parabolic 
trend  of  wave  period  that  was  observed  in  the  data  for  slopes  of  1  on  2.5  and 
1  on  5,  as  discussed  by  Ahrens  and  McCartney  (1975),  while  the  small-scale  Nz's 
show  a  more  linear  trend. 


Profile  Shapes. 


Another  way  to  evaluate  scale  effects,  particularly  at  high  damage  levels, 
is  to  compare  the  shapes  of  the  damaged  surface  profiles  for  similar  wave 
conditions.  This  comparison  requires  the  use  of  data  where  the  dimensionless 
damage  is  about  the  same  at  both  small  scale  and  prototype.  Figure  8  shows  a 
profile  comparison  for  tests  with  short  period  waves  (d/gT2  ■  0.0264)  and  gives 
the  dimensionless  wave  height  and  damage,  respectively,  for  both  tests  in  the 
legend.  The  profiles  were  made  comparable  by  increasing  the  small-scale  test 
dimensions  by  a  factor  of  10.  The  shapes  of  the  two  profiles  were  similar,  and 
the  causative  dimensionless  wave  heights  were  about  the  same.  Figure  9  is 
similar  to  Figure  8  except  the  profile  comparison  was  for  tests  with  long 
period  waves  (d/gT2  ■  0.0037),  and  the  causative  dimensionless  wave  height  was 
approximately  18  percent  smaller  in  the  small-scale  tests. 
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The  correlation  coefficients  of  the  damage  profiles  provided  a  convenient 
means  of  comparing  the  damage  profiles  of  the  model  and  the  prototype  by  cal¬ 
culating  the  coefficients  for  model  and  prototype  tests  with  approximately  the 
same  relative  damage.  The  profiles  tests  were  matched  by  the  location  of  the 
Stillwater  level  on  the  reference  surveys,  and  the  averaged  differences  from 
the  reference  surveys  for  the  model  and  prototype  tests  were  paired  to  calcu¬ 
late  the  correlation  coefficient.  Table  2  provides  a  tabulation  of  correlation 
coefficients  for  a  number  of  model-prototype  profile  pairs,  including  those 
shown  in  Figures  8  and  9.  In  general,  the  closer  the  relative  depths  of  the 
model  and  prototype  tests  were  the  more  similar  the  profile  shapes. 

3.  Runup . 

In  evaluating  scale  effects  between  small-scale  and  large-scale  tests  the 
wave  runup  was  also  compared.  Runup  was  visually  defined  as  the  average  point 
of  maximum  wave  uprush  on  the  riprap  surface  near  the  center  of  the  wave  tank. 
The  elevation  of  this  point  was  then  measured  using  the  survey  apparatus.  The 
runup  data  in  Table  3  indicate  that  relative  runup,  R/H,  can  be  considered 
constant  for  a  fixed  value  of  d/gT2  between  0.0264  and  0.0036.  For  d/gT2 
of  0.0595  or  0.0589,  the  ratio  of  the  relative  runup  to  the  surf  parameter,  C, 
is  almost  constant,  where  the  surf  parameter,  £,  is  defined  as 

t  «  (H/Lo)’1/2  tanB, 

L0  is  the  deepwater  wavelength,  and  Q  the  angle  between  the  embankment  and 
the  horizontal.  The  runup  invariants,  R/II  or  (R/H)//;,  as  tabulated  in 
Table  3  and  shown  in  Figure  6,  indicate  that  the  runup  in  the  small-scale  tests 
was  approximately  20  percent  greater  than  predicted  from  the  prototype  tests. 

The  small-scale  test  results  at  the  zero-damage  level  give  more  conservative 
estimates  of  the  stability  number  and  the  runup.  The  zero-damage  level  stabil¬ 
ity  numbers  are  lower  and  the  runup  is  higher  in  the  small-scale  tests  than 
predicted  by  the  prototype  tests.  The  higher  runup  is  probably  due  to  the 
reduced  penetration  of  the  wave  uprush  In  the  small-scale  tests  as  compared  to 
the  prototype  tests.  The  stone  size  used  in  the  filter  layer  was  modeled 
geometrically  but  should  be  somewhat  larger,  according  to  Keulegan  (1973),  to 
obtain  proper  flow  similitude. 

4 .  Flow  Regime. 

In  this  study  the  small-scale  tests  replicated  the  prototype  tests  at  a 
1:10  Froude  scale  which,  assuming  no  scale  effects,  required  rough  turbulent 
flow  in  both  the  small-scale  and  prototype-scale  tests.  When  using  the  Froule 
scale  model  the  model  and  prototype  will  be  dynamically  similar  with  respect 
to  inertial  and  viscous  forces  because  viscous  forces  can  be  assumed  insignifi¬ 
cant.  The  existence  or  nonexistence  of  rough  turbulent  flow  is  determined  from 
the  criterion  established  by  Jonsson  (1966).  Using  definitions  similar  to 
Madsen  and  White  (1976),  who  applied  Jonsson's  criterion  to  rubble-mound 
structures,  the  Reynolds  number,  RE,  is  given  by 

F2  (1  +  cot20)  (|— ) 

r£ - ±—  (5) 
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Table  2.  Correlation  coefficients  for 
model  and  prototype  profiles. 


Model 

Prototype 

Test 

SPL  25 

SPL  19 

SPL  22 

designation 

d/gT2 

0.0264 

0.0144 

0.0037 

SET-2 

0.0264 

0.847 

0.927 

0.796 

SET-1 

0.0144 

0.848 

0.899 

0.752 

SET-13 

0.0037 

0.607 

0.702 

0.946 

R/K 

Mean 

Std.  dev. 

0.682 

0.760 


0.024 

0.051 


0.0064  1.444 
0.0064  1.649 

0.0037  1.568 
0.0037  1.948 


0.040 

0.033 

0.021 


0.038 

0.097 


and  the  roughness  term,  aim/k, 


by 


aim  Rz  d  +  cot^)1/2 


(6) 


where  Rz  is  the  wave  runup  associated  with  the  zero-damage  wave  height,  Hz; 
T  the  wave  period;  and  v  the  kinematic  viscosity.  The  values  of  runup  used 
to  compute  the  Reynolds  number  and  roughness  are  tabulated  by  test  in  Table  4 
and  represent  the  estimated  runup  which  would  be  caused  by  the  zero-damage 
wave  height.  The  estimates  of  Rz  are  calculated  using  the  values  of  Hz  in 
Table  1  with  the  runup  invariants  tabulated  in  Table  3.  Since  the  Reynolds 
number  defined  in  equation  (5)  uses  wave  runup,  the  calculations  of  the  flow 
regime  refer  to  surface  conditions,  not  conditions  in  the  filter  layer.  In 
Figure  10  the  Reynolds  number  and  roughness  values  tabulated  in  Table  4  are 
shown  with  the  flow  regime  boundaries  as  updated  by  Jonsson  (1978) .  The 
figure  shows  that  both  the  small-scale  and  prototype  tests  are  in  the  rough 
turbulent  flow  regime. 


V.  COMPARISON  WITH  OTHER  SOURCES  OF  DATA 

Scale  effects  at  the  zero-damage  level  were  compared  with  the  scale  effects 
test  results  of  Dai  and  Kamel  (1969)  and  Thomsen,  Wohlt,  and  Harrison  (1972)  in 
Figure  11;  Dai  and  Kamel  used  rough  quarrystone  in  their  rubble-mound  stability 
tests  and  Thomsen,  Wohlt,  and  Harrison  used  dumped  Kimmswick  limestone  in  their 
riprap  stability  tests.  The  comparison  was  made  by  dividing  the  individual 
value  of  Nz  by  the  average  prototype  value  of  Nz  for  the  tests  having  the 
same  relative  depth  (designated  NZp^,  as  tabulated  in  Table  1  for  this  study, 
to  form  the  scale  effects  factor  Nz/Nzp.  Figure  11  shows  the  scale  effects 
factor  plotted  versus  a  Reynolds  number,  Rfj,  which  is  given  by 


where  Hz  is  the  zero-damage  wave  height,  and  v  the  kinematic  viscosity  of 
water  (assumed  to  be  1.1306  x  10-6  square  meters  per  second  corresponding  to 
a  water  temperature  of  15.6°  Celsius).  At  the  zero-damage  level  the  stability 
numbers  were  approximately  20  percent  lower  for  the  small-scale  test  than  for 
the  prototype  test,  as  shown  in  Figure  11.  The  figure  also  shows  that  at  the 
zero-damage  level  the  scale  effects  observed  in  this  study  were  somewhat  less 
severe  than  those  observed  by  Thomsen,  Wohlt,  and  Harrison  and  comparable  to 
those  observed  by  Dai  and  Kamel. 

The  small-scale  test  results  of  this  study  were  also  compared  with  a  1975 
study  conducted  by  the  Hydraulic  Research  Station  (HRS),  Wallingford,  England, 
for  the  Construction  Industry  Research  and  Information  Association  (CIRIA) 
(Hydraulic  Research  Station,  1975).  The  HRS  study  on  riprap  stability  under 
irregular  wave  attack  was  conducted  at  small  scale.  Table  5  tabulates  the  HRS 
tests  with  a  1  on  4  or  1  on  3  slope  which  were  used  for  comparison.  These  slopes 
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Table  4.  Flow  regime  computations-. 


Test 

d  /  gT  “ 

T 

'■'5C 

. 

V 

k 

Rz 

aim/k 

rF. 

designation 

1  X  10~5 

1  x  10'- 

(eq_.  6) 

(e<q.  5) 

SPL-17 

0.0595 

2.8 

75 

0.97 

0.763 

3.48 

16.60 

371 

SET-22 

0.0588 

0.89 

0.075 

1.06 

0.0763 

0.313 

14.94 

8.65 

SET-24 

0.0588 

0.89 

0.075 

1.02 

0.0763 

0.335 

15.98 

10.29 

SPL-1S 

0.0264 

4.2 

75 

0.92 

0.763 

3.21 

15.32 

222 

SPL-25 

0.0264 

4.2 

75 

1.20 

0.763 

2.86 

13.64 

135 

SET- 2 

0.0264 

1.33 

0.075 

0.99 

0.0763 

0.300 

14.33 

5.71 

SET-9 

0.0264 

1.33 

0.075 

1.003 

0.0763 

0.264 

12.59 

4.36 

SET-19 

0.0264 

1.33 

0.075 

1.00- 

0.0763 

0.234 

11.15 

3.42 

SPL-19 

0.0144 

5.7 

75 

1.08 

0.763 

2.93 

14.00 

116 

SPL-23 

0.0144 

5.7 

75 

1.50 

0. 763 

2.67 

12.75 

63 

SET-1 

0.0144 

1.80 

0.075 

1.00- 

0.0763 

0.264 

12.58 

3.21 

SET- 3 

0.0144 

1.80 

0.075 

1.01 

0.0763 

0.251 

11.93 

2.89 

SET- 10 

0.0144 

1.80 

0.075 

1.003 

0.0763 

0.256 

13.62 

3.77 

SET-20 

0.0144 

1.80 

0.075 

1.13 

0.0763 

0.291 

13.82 

3.43 

SPL-20 

0.00646 

8.5 

75 

1.75 

0.763 

3.86 

18.41 

83 

SPL-24 

0.00646 

8.5 

75  • 

1.48 

0.763 

3.80 

18.13 

96 

SET- 6 

0.00646 

2.69 

0.075 

1.06 

0.0763 

0.386 

18.41 

4.35 

SET- 11 

0.00646 

2.69 

0.075 

1.003 

0.0763 

0.388 

18.49 

4.65 

SET-21 

0.00646 

2.69 

0.075 

1.12 

0.0763 

0.388 

18.49 

4.15 

SPL-22 

0.00365 

11.3 

75 

1.59 

0.763 

4.87 

23.26 

110 

SET-12 

0.00366 

3.57 

0.075 

1.003 

0.0763 

0.425 

20.26 

4.21 

SET-13 

0.00366 

3.57 

0.075 

1.003 

0.0763 

0.419 

19.98 

4.09 

SET-23 

0.00366 

3.57 

0.075 

1.03 

0.0763 

0.506 

24.16 

L_!\21.__.| 

3A11  test  cot@  =  3.5. 


2k  =  (W50/wr>  . 

3As-sumed  value  of  kinematic  viscosity. 
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Table  5.  HRS  data. 


Tests 

Slope 

d 

(cm) 

■ 

■ 

1 

b  ‘ 

h2  2 

(cm) 

Rz  3 

(cm) 

re  4 
1  *  104 

49-61 5 

1:3 

61.0 

1.64 

1.1040 

2.7063 

3.08 

4.23 

72.9 

108-1126 

1:3 

30.5 

1.64 

1.4184 

5.0614 

3.95 

5.98 

209.2 

33-46 5 

1:3 

61.0 

2.46 

1.2455 

3.2222 

5.21 

8.49 

29.4 

113-1167 

1:3 

46.1 

2.46 

1.4122 

4.7379 

5.91 

9.88 

460.9 

62-67 5 

1:3 

61.0 

3.28 

1.3621 

3.5014 

7.50 

13.26 

71.8 

117-127  5 

1:4 

61.0 

2.0 

1.64 

1.6297 

4.2240 

4.54 

5.41 

203.5 

166-169® 

1:4 

30.5 

2.0 

1.64 

1.7601 

4.0394 

4.91 

5.87 

342.8 

128-1405 

1:4 

61.0 

IB9 

2.46 

1.6869 

5.4265 

7.05 

8.55 

507.7 

170-1738 

1:4 

61.0 

m 

3.28 

1.6610 

5.4753 

9.26 

11.31 

888.0 

«  D 

•  An  equation  of  Che  form  D'  -  aNs  was  fitted  to  the  data  . 

2 

Hz  predicted  wave  height  at  D'  =■  2.00. 

3Rz  predicted  runup  for  Hz  . 


Re  *  - ^ -  ,  where  v  *  9.3  *  10  3  square  centimeters 

per  second  and  d/gT2  =  0.023  is  assumed. 

5Tp  ranged  from  1.66  to  1.15  seconds,  d/gT2  ranged  from  0.023  to  0.047. 
6Tp  equaled  1.154  seconds,  d/gT2  equaled  0.023. 

7Tp  equaled  1.431,  d/gT2  equaled  0.023. 

3Tp  equaled  1.655,  d/gT2  equaled  0.023. 
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were  chosen  because  they  bracketed  the  1  on  3.5  slope  used  in  the  small-scale 
tests  of  this  study.  Only  the  small-scale  tests  with  a  relative  depth  of 
0.0144  and  0.0264  were  used.  The  relative  depth  of  0.0144  gave  the  lowest 
stability  numbers  and  0.0264  was  in  the  range  of  relative  depths  tested  in  the 
HRS  study.  Figure  12  shows  the  stability  numbers,  Nz,  versus  the  Reynolds 
number,  Rg,  which  is  described  by  equation  (5).  From  Figure  12  it  can  be 
hypothesized  that  riprap  stability  under  irregular  wave  attack  may  be  equiva¬ 
lent  to  the  monochromatic  tests  with  the  relative  depths  that  yield  the  lowest 
stability  numbers. 

VI.  RESULTS  AND  CONCLUSIONS 

The  results  of  this  study  showed  a  reduction  of  about  20  percent  in  the 
zero-damage  stability  numbers  for  a  1:10  (model prototype)  Froude  scale  model 
from  the  expected  prototype  values.  The  reduction  of  stability  in  the  model 
appears  to  be  related  to  the  lack  of  penetration  of  the  wave  runup  into  the 
filter  layer  and  the  improper  modeling  of  the  flow  regime  within  the  filter 
layer.  Figures  5  and  7  show  that  the  difference  between  the  small-scale  tests 
and  the  prototype  tests  decreased  as  the  damage  level  increased  indicating 
that  the  scale  effects  decrease.  The  data  points  in  Figure  7  indicate  a 
convergence  of  the  damage  trends,  whereas  the  equation  D'  =  aN|,  as  shown  in 
Figure  5,  shows  a  crossing  of  the  damage  trends.  The  convergence  of  damage 
trends  seems  reasonable  because  higher  Reynolds  numbers  developed  at  higher 
damage  level  and  the  viscous  forces  became  less  significant.  The  crossing  of 
the  damage  trends  seems  unreasonable  and  could  be  caused  by  the  method  used  in 
deriving  the  equation  D'  =  aN^.  The  data  points  at  the  lower  damage  levels 
could  be  exerting  more  influence  on  the  equation  than  is  justified. 

The  following  conclusions  were  reached: 

1.  The  tests  at  a  1:10  (model : prototype)  Froude  scale  yield  zero-damage 
stability  numbers  about  20  percent  lower  than  the  prototype  tests.  This 
indicates  that  scale  effects  in  this  study  were  less  severe  than  those  found  by 
Thomsen,  Wohlt,  and  Harrison  (1972). 

2.  Scale  effects  were  less  severe  at  high  levels  of  damage  than  at  the 
zero-damage  level. 

3.  The  runup  at  the  zero  level  was  about  20  percent  higher  in  the  small- 
scale  tests  than  predicted  by  the  prototype  test. 

4.  The  shapes  of  the  damage  profile  for  the  small-scale  and  the  prototype 
tests  having  the  same  relative  depth  were  very  similar. 

5.  At  the  zero-damage  level,  wave  period  had  less  influence  in  the  small- 
scale  tests  than  in  the  prototype  tests. 
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HRS  I  on  3  Slope 
HRS  I  on  4  Slope 
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